Late-time Kerr tails revisited 
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The decay rate of late time tails in the Kerr spacetime have been the cause of numerous conflicting 
results, both analytical and numerical. In particular, there is much disagreement on whether the 
decay rate of an initially pure multipole moment I is according to t~^ 2e+3 \ where I is the least 
multipole moment whose excitation is not disallowed, or whether the decay rate is according to t~ n , 
where n = n(£). We do careful 2+1D numerical simulations, and explain the various results. In 
particular, we show that pure multipole outgoing initial data in either Boyer-Lindquist on ingoing 
Kerr coordinates on the corresponding slices lead to the same late time tail behavior. We also show 
that similar initial data specified in terms of the Poisson spherical coordinates lead to the simpler 
£-(2i+3) | a ^ e tj me £ a ji -y^g generalize the rule n — n(l) to subdominant modes, and also study 
the behavior of non-axisymmetric initial data. We discuss some of the causes for possible errors in 
2+1D simulations, demonstrate that our simulations are free of those errors, and argue that some 
conflicting past results may be attributed to them. 



PACS numbers: 04.70.Bw, 04.25.Nx, 04.30.Nk 

I. INTRODUCTION AND SUMMARY 

The late-time tails of black holes have been studied 
in much detail since Price's seminal work pj. The for- 
mulation of the problem is a straightforward one: place 
an observer in a circular orbit around a black hole, and 
have her measure at late times a generic perturbation 
field, that had compact support at some initial time. It 
is generally accepted that the observer measures the late- 
time perturbation field to drop off as an inverse power law 
of time, specifically as t~ n . It is the value of n that has 
been controversial in the literature, with some conflicting 
results reported. 

In the case of a Schwarzschild black hole, n — 21 + 3, 
where £ is the multipole moment of the initial perturba- 
tion field. Namely, if the initial (compactly supported) 
perturbation field has the angular dependence of Y™ , the 
angular dependence remains unchanged ("spherical har- 
monics are eigenvectors of the Laplacian operator" ) , and 
the decay rate of the field is governed by the £ value of 
the initial perturbation. These results remain unchanged 
also for a Rcissner-Nordstrom black hole (includ mg 
the extremal case Q), that like the Schwarzschild black 
hole is spherically symmetric. Notably a generic pertur- 
bation is a linear superposition of £-modes, so that the 
tail is dominated by the slowest damped mode, or by the 
lowest value of £. 

For rotating black holes conflicting values of n have 
been reported. Obviously, the observer measures a well 
defined decay rate. Of the various conflicting values in 
the literature some are wrong (i.e., they report a value 
inconsistent with the outcome of the initial value prob- 
lem), while other results are correct (in the same sense, 
of reporting the right value for the initial value problem) . 
As it turns out, different initial value formulations may 



lead to different decay rate for the tails. While none of 
these is "more correct" sensu stricto than others, it is im- 
portant to understand the relationship of different initial 
value sets, so that the understanding may lead to some 
insights. 

One of the correct answers to the question of the value 
of n was described, e.g., in 12]: the late-time field de- 
cays according to the same rule as in the spherically- 
symmetric case, namely according to i~( M+3 ), where 
£ = to if I — m is even, and £ = m + 1 if £ — to is odd, in- 
dependently of the initial value of £. Below, we dub this 
tail as "Tail A." The reason for this decay rate is that 
spherical harmonic modes do not evolve independently 
on the background of a rotating black hole ("spherical 
harmonics are not eigenvectors of the Laplacian opera- 
tor"). Starting with an initial value for £, infinitely many 
other £ modes are excited, respecting only the dynamical 
constraints (£ > s, where s is the spin of the field, and 
£ > | to |) and the equatorial symmetry of the initial data 
(such that even and odd modes do not mix). 

Another correct answer to the question was discussed 
recently in 0] and in According to this result, n — 
2£ + 3 if £ - to < 2, n = £ + m+ 1 if I - to > 2 is 
even, and n = £ + m + 2 if £ — to > 2 is odd. This 
result is dubbed below as "Tail B." Evidently, this result 
is inconsistent with that of Tail A. The Tail B behavior 
was first predicted in While either the Tail A or 
the Tail B behavior is a solution to some initial value 
problem, these problems are different enough to lead to 
markedly different late-time evolutions. 

In this paper we present some explanation for these 
different results, that resolves the apparent discrepancy 
at least for some cases. This paper's motivation is to ex- 
plain some of the historical discrepancies, in what we be- 
lieve is an illuminating manner. Specifically, finding the 
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reasons for the controversy in the literature teaches us 
important lessons. In particular, we consider two results, 
conflicting both with each other and with the results of 
this paper: the numerical results of [H, and the analyti- 
cal model of [|| . More precisely, we consider the simplest 
case for which there is a disagreement in the literature, 
specifically the case of a massless scalar field (s = 0) per- 
turbation of a Kerr black hole, with initial I = 4 and 
to = 0. For this case, the late time tail drops as t~ 3 ac- 
cording to Tail A, because the monopole (£ = 0) mode of 
the scalar field can and is excited, and this mode's drop 
off rate obeys the i~( 2£ + 3 ) rule. According to Tail B, the 
late-time decay rate is i~ 5 . The (numerical) results of 
H suggest a drop off rate of £~ 5 - 5 that does not appear 
to converge to the decay rate of Tail B. The Poisson an- 
alytical model (of a globally weakly curved spacetime) 
predicts a drop off rate of t~ 5 if carried over to the Kerr 
case Q, which may suggest a Tail B behavior. How- 
ever, as we show below, the correct behavior of such tails 
should be that of Tail A. 

We discuss in detail the tails developing from initial 
data specified on ingoing Kerr silces oflj]. It was sug- 
gested in Q that the Tail A results of [j] are due to the 
slicing condition. We show that this is not the case. In 
fact, one can find initial condition sets with the same 
slicing condition as in J| (slicing corresponding to the 
ingoing Kerr coordinates) that result in either Tail A or 
Tail B decay rates, in just the same way as with Boyer- 
Lindquist slicing. In fact, although the Boyer-Lindquist 
and ingoing Kerr time slices do not coincide, we find 
that similar initial data lead to the same tail behav- 
ior in cither. We therefore suggest that different slicing 
conditions may fall into equivalency classes, and Boyer- 
Lindquist and ingoing Kerr slices belong to the same 
class. We propose that this is the case because of the 
properties of the transformation from Boyer-Lindquist 
and ingoing Kerr coordinates. 

In the case of Q, we attribute the (numerically sta- 
ble!) result of a drop off rate of t~ 5 - 5 to an intermedi- 
ate tail being misidentified as an asymptotic tail. As we 
show, the traditional representation of the tail in a log- 
log plot of the field as a function of time is susceptible 
to lead to an "optical illusion," that may support this 
misidentification. Specifically, taking initial conditions 
that lead to a Tail A decay rate, we show what may be 
interpreted as a tail falling off as i~ 5 - 5 , similar to what is 
found in Q, except that if the integration is carried on 
further in time the intermediate behavior is taken over 
by the true asymptotic behavior, or a tail falling off as 
i -3 . A similar effect can also be found for the Tail B 
case. The reason Q mistook an intermediate tail for an 
asymptotic one is that [f| chose parameters that make it 
harder to distinguish the two using the crude tools used 
in Specifically, we show that the farther out the ini- 
tial pulse, the longer the intermediate tail. Reference Q 
chose the initial pulse very far out, which resulted in a 
protracted intermediate tail. While this choice made the 
problem severe enough, it was even further exacerbated 



by making the initial pulse outgoing and by taking the 
Kerr parameter a/M = 0.9999. Either choice makes the 
beginning of the tail regime postponed. As the integra- 
tion time in [5[ is rather short, we conclude that the tail 
analyzed in [5[ is an intermediate one, and that if the in- 
tegration time were longer, [f| too would have found a t~ 3 
tail (for Tail A initial data) or a t~ 5 tail (for Tail B initial 
data) . (In view of our argument below, we comment that 
the code of [|[ was fourth order angularly, which we find 
to be satisfactory for the relatively short runs done by 
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How can one then avoid misinterpreting intermediate 
tails as asymptotic ones? We provide a criterion, that we 
believe can be used for the identification of the asymp- 
totic tail. By our criterion, a tail is asymptotic if the local 
power index Q is linear in i -1 both locally and globally, 
in a sense that we make precise. 

In the case of @, we explicitly follow the proposal of 
Poisson, and evolve initial data that are a pure spher- 
ical harmonic mode with respect to "spherical" coor- 
dinates, not the Boyer-Lindquist coordinates, that are 
"spheroidal" in the language of Q . According to [§:] , such 
initial data should lead to a Tail B behavior when applied 
to a Kerr background. We show that this is not the case, 
and Tail A behavior is found. Specifically, starting with 
I = 4 and m — 0, Poisson predicts the asymptotic tail 
to drop off as i -5 . In [6(, such time evolution is demon- 
strated in a globally weak-curvature spacetime (that has 
no symmetries), and it is conjectured that such behavior 
is carried over to the case of a Kerr black hole. We show 
that this conjecture is not realized: instead, the evolution 
is according to Tail A, and the drop off rate is according 
to t~ 3 . 

Very recently, Gleiser, Price, and Pullin (henceforth 
GPP) have proposed a theoretical model that expands 
the wave equation in (even) powers of the black hole's 
spin rate a [9j . The advantage of the GPP approach is 
that the zeroth order wave operator is identical to the 
spherically-symmetric operator, so that it can be solved 
numerically using a simple 1+1D code. Each iteration 
in the GPP approach gives the field to the next order in 
a 2 , using the field of the previous iteration as the source 
term of a spherically symmetric wave operator. GPP 
argue that their approach is superior to the 2+1D ap- 
proach, because it frees the numerics from the subtleties 
of the 2+1D simulations, and because the 1+1D numerics 
are inherently more reliable. We agree with these claims. 
The GPP argument implies that delicate cancellations 
exactly cancel not only the leading tail term (in an ex- 
pansion in t~ l ), but also the leading subdominant mode 
[llT | (for an initial £ = 4; for other initial choices more 
subdominant modes may need to be cancelled). GPP 
further argue that turning off the coupling term at very 
late times destroys the exact cancellation, so that the tail 
returns to the "naive" decay rate. It is quite remarkable 
that even if this is done at extremely late times, so that 
the hexadecapole field is much smaller than the monopole 
field, it still does not lose its ability to exactly cancel out 
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its two leading terms (in t 1 ). 

Notably, GPP imply that different slicing conditions 
may be responsible for the apparently conflicting results 
of numerical simulations. Indeed, in [4J the slicing condi- 
tion was that associated with ingoing Kerr coordinates, 
in [l2j the slicing was that of the Kerr-Schild coordi- 
nates, and in [j| Boyer-Linquist slicing was used. GPP 
argue that the translation of an initially pure I = 4 mode, 
say, in some other slicing condition (e.g., ingoing Kerr 
or Kerr-Schild) into Boyer-Lindquist slicing results in 
a mixed state of modes for the Boyer-Lindquist slicing, 
and the existing lower £ (specifically, the quadrupole and 
monopole) modes on the initial time slice are the ones 
responsible for the tail structure. According to GPP, 
among all slicing conditions the Boyer-Lindquist one is 
special, as only in this slicing an initial pure I — 4 mode 
leads to a tail of t~ 5 , and in all other slicing conditions an 
initially pure £ = 4 mode leads to a t~ 3 tail. We disagree 
with this viewpoint, as we find that also for ingoing Kerr 
slicing the tail decays like t~ 5 . 

More recently, Tiglio, Kidder, and Teukolsky (TKT) 
[j~0| performed careful 2+1D simulations of the axisym- 
metric Kerr tail problem using a pseudospectral code. 
The pseudospectral approach has an advantage over the 
finite differencing approach, in that the convergence is 
exponential in the former, and only a power law in the 
latter. Consequently, in the pseudospectral approach the 
noise evolution problem that permeates finite differenc- 
ing computations is nearly nonexistent. TKT studied two 
slicing cases: Kerr-Schild slicing and Boyer-Lindquist 
slicing. TKK found the Tail A behavior in the former 
(consistently with the conclusions of [l2j that were ob- 
tained in 3+ ID) and the Tail B behavior in the latter 
(consistently with GPP). TKT correctly point out the 
main reason for the difference between Tail A and Tail B 
behaviors that are found in Boyer-Lindquist initial data 
and on Boyer-Lindquist slicing, specifically that the zero 
"field momentum" of the former is in fact a mixed state 
of pure multipole modes. We consider this question in 
great detail, and show that the momentum condition is 
a very special combination of states, as described below. 

We discuss in detail the noise evolution problem that 
may cause a Tail B masquerade as a Tail A. That is, 
numerical noise evolves on the grid the same way that 
the physical signal does, so that the Tail A behavior of 
the noise evolution overwhelms the Tail B behavior of 
the physical field. We then show that with care this 
phenomenon can be easily identified. 

We next discuss how initial data with the same slicing 
condition could lead to either Tail A or Tail B behavior. 
As shown first by GPP and then by TKT, initial data 
that are a pure multipole of the angular Boyer-Lindquist 
coordinates in Boyer-Lindquist slicing lead to Tail B be- 
havior. When one sets the initial conditions by requiring 
that the field's "momentum" n (see below for definitions) 
vanishes, the tail's behavior is that of Tail A. Indeed, the 
condition on the momentum implies that the field's mul- 
tipole structure is changed, as the field's velocity now is 



proportional to the field multiplied by a function of the 
angular coordinate 9. Restricting attention to low black 
hole spin values (a/M < 1), this function is proportional 
to a 2 sin 2 6. Therefore, one would expect that in addition 
to an initial hexadecupole at O(a ) one would have now 
also all other (even) multipoles at 0(a 2 ). An existing 
monopole on the initial slice at 0(a 2 ) would then lead to 
a tail that has the same decay rate as a pure multipole 
that behaves like Tail A, which is at 0(a 2 ). However, 
we show that the tail is at 0(a 4 ), not 0(a 2 ). The ex- 
planation is that the leading contribution (in a 2 ) does 
not contribute to the monopole, so that on the initial 
slice the monopole is at 0(a 4 ), not at 0(a 2 ). Evolution 
of this mode then can dominate at late times, showing 
a tail that is both at 0(a 4 ) and has the Tail A behav- 
ior. Nevertheless, this is not the only contribution to this 
tail. The initial data contains also a quadrupole at 0(a 2 ), 
which also evolves a 0(a A ) monopole because of the mode 
excitation effect. However, as in this case the quadrupole 
has an £ value that differs from the monopole's by only 
2, the behavior is again consistent with that of Tail A. 
Therefore, the Tail A behavior has three sources: a pre- 
existing monopole at (3(a 4 ), a pre-existing quadrupole 
at 0(a 2 ), and a pre-existing hexadecapole at 0(a°), all 
three behaving consistently with Tail A. 

We further extend the numerical study to non- 
axisymmetric modes, and find that the Tail A and Tail 
B behaviors occur for such modes consistently with their 
behavior in the axisymmetric case, i.e., they obey the 
aforementioned decay rate rules. We then extend the 
argument for Tail B behavior to subdominant modes. 
Specifically, we argue that starting with an initial mode 
£, the decay rate of an (allowed) excited mode £' is given 
by n = I + 1' + 1 if I - I' > 2, and n = 21' + 3 otherwise. 
This rule includes that of Q as a particular case, when 
£' = m^m + 1. In particular, we consider the axisym- 
metric case of £ = 6 and £' — 2 (which is a subdominant 
mode, as also the monopole is excited). The quadrupole 
mode £' = 2 then indeed decays with n = 9 as predicted 
by our proposal. 

We are hopeful that this Paper may put some order in 
a confused state of affairs, in which the literature includes 
side by side correct results and incorrect ones, without a 
proper explanation of the incorrect ones, and also appar- 
ently conflicting correct results without a proper expla- 
nation of the sources for the apparent conflict. 

The organization of this paper is as follows. In Sec- 
tion [TT] we discuss the reliability of the 2+ ID approach 
in calculating tails, and how carelessness can result in an 
incorrect value for the field. In Section HTT1 we discuss the 
decomposition of the monopole projection of the tail in 
powers of a 2 , and in Section [V] we discuss one possible 
pitfall that can cause a correct field be misinterpreted. 
In Section IVT1 we study tails in ingoing Kerr coordinates, 
and in Section [VlT] we analyze a pure mode in "spherical" 
coordinates [6|], and show that it leads to the same tail 
decay rate as in "spheroidal" coordinates. 
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II. RELIABILITY OF THE 2+1D CODE 

The solution of the Teukolsky equation in 2+1D is del- 
icate and sensitive, and can lead to qualitatively wrong 
results if proper care is not taken. In Section Ivl below we 
show that even when the numerically computed field is 
the correct one, careless interpretations may lead to in- 
correct conclusions. In this section we first demonstrate 
how carelessness may lead to qualitatively wrong behav- 
ior of the solution, and then we demonstrate the reliabil- 
ity of the approach when proper care is taken. 



A. Second order angular differentiation 

The standard of scientific computation in the finite dif- 
ference approach has been second order algorithms for 
many years. Indeed, second order convergence typically 
implies the solution is both consistent and stable to a 
practical level, so that the solution is accurate enough 
without having to reduce the typical grid spacing to ex- 
tremely low levels that would make the code too costly 
computationally. Normal application of the Dahlquist 
theory (extended to PDEs) then means that if the dis- 
crete formula approximating the differential equation has 
locally a positive order of accuracy (convergence), and 
that if numerical errors do not grow unboundedly (sta- 
bility), then the solution is convergent, in the sense that 
as the grid spacing goes to zero, the solution would ap- 
proach the correct one. As we show below, the tail prob- 
lem involves the Dahlquist theory with a twist. 

Early codes to calculate the late-time Kerr tail were 
therefore designed to be second-order in all coordinates, 
i.e., temporally, radially, and with respect to the angu- 
lar coordinate 9 [HI]. (In the 2+ ID approach the field 
is decomposed into ip modes, so that one does not have 
numerical derivatives with respect to (p.) While this ap- 
proach seems very natural, it can mislead into qualita- 
tively wrong solutions. Consider first a Schwarzschild 
black hole, excited by an azimuthally symmetric (m = 0) 
field with t = 4. The decay rate of the Price tail is in- 
controvertible in this case: as the background is spheri- 
cally symmetric, spherical harmonic modes evolve inde- 
pendently, and the decay rate is therefore t~ xx . Indeed, 
I+1D codes readily verify this result. 

The numerical solution with a second-order 2+ ID nu- 
merical code may lead however to the erroneous conclu- 
sion that the decay rate is t~ 3 : in Fig. [T] we show the 
late-time field as a function of time for three choices of 
the angular grid spacing, and for fixed radial and tem- 
poral grid spacings. As can be readily inferred from the 
upper panel in Fig.[TJ the decay rate is t~ 3 , which is defi- 
nitely wrong. The origin of this wrong result is indicated 
in the same figure: the field values themselves reduce 
considerably in amplitude under refinement of the grid 
(almost as the second power of the ratio of the grid pa- 
rameters). The signal seen is therefore not the physical 
field, but an evolved numerical noise ("leakage of multi- 
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FIG. 1: Convergence tests for Schwarzschild initial data with 
the second-order code. Upper panel and lower panel on the 
left: the full field (evaluated at 6 — n/2) for three different 
angular resolutions (A cost* = 1/72 (solid), 1/48 (dashed), 
and 1/32 (dotted)) as a function of time. The lower panel on 
the right is the convergence order as a function of the time 
at r*/M = 25. The initial data are for a Schwarzschild black 
hole (a/M — 0), and are a gaussian centered at r* = 25M 
with width AM. The radial resolution is Ar = Af/32, and 
the temporal resolution is chosen to be At = Ar/2, which 
satisfies the Courant condition. The initial data are a pure 
multipole with t — 4. 



poles" ) . Our interpretation is that while the physical sig- 
nal indeed drops like t , it is swamped by the evolved 
numerical noise, which has a monopole component. The 
latter drops like i -3 , so that any very long time evolu- 
tion must be dominated by numerical noise. Refinement 
of the grid parameter can postpone the numerical-noise 
dominated regime, but cannot utterly eliminate it: re- 
finement of the grid parameter reduces the level of the 
numerical noise, so that if the grid is fine enough, the 
noise level at a particular value of the time reduces to 
below the level of the physical field. However, as the for- 
mer drops more slowly than the latter, at very late times 
the field is dominated by numerical noise (see Fig. [5] be- 
low). 

Here enters the twist on the Dahlquist theory: the nu- 
merical noise does not grow unboundedly. In fact, it even 
drops quadratically with the grid spacing, and drops with 
time. It does so, notwithstanding, slower than the phys- 
ical signal, so that the signal may be dominated by nu- 
merical noise even though the regular Dahlquist criteria 
are apparently satisfied. 



B. Sixth order angular differentiation 

While a second order code may in principle be sufficient 
to see physical tails (for very fine grids), an alternative 



5 



approach is to increase the order of the angular differ- 
entiation. We first show how the same initial data and 
grid resolutions lead to qualitatively different results in 
Schwarzschild, and then we discuss Kerr. 

Figure [2] shows the sixth order convergence of the field, 
for the same parameters that were shown in Fig. [T] with 
the second order code. The field is not made of evolved 
noise ("leakage of multipolcs" ) , as it clearly converges 
to a non-zero value under grid refinement. Notably, the 
field's amplitude even increases slightly under grid re- 
finement. The hallmark of noise signals is that their am- 
plitude decreases under grid refinement, as is indeed ev- 
ident in Fig. [TJ Most importantly, the Price tail drops 
as t -11 , as it should. We emphasize that these conclu- 
sions are true up to some maximal time value. For longer 
times the residual monopolc noise, that drops slower than 
the physical field, will inevitably dominate. Simple ex- 
trapolations reveal that the signal noise will dominate 
over the physical signal for t > 4, 380Af for our choice 
of angular resolution. Longer evolutions than that will 
be dominated by noise evolution. As noted above, it is 
not difficult to separate the physical signal regime from 
the noise signal regime by doing global convergence tests. 
It is therefore crucial for any tail simulation that conver- 
gence tests are done throughout the entire computational 
domain, and not limited just to short integration times 
(despite the temptation to save computer time — e.g., 
see (|). 
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FIG. 2: Same initial data and resolutions as in Fig. [T]for the 
sixth order code. The upper panel shows the L2 norm for 
the full field (integrated over a sphere at r* = 25M) in three 
resolutions. The three curves are nearly indistinguishable on 
this plot's scale. The upper plot also shows the monopole 
projection of the field, which is pure evolved numerical noise 
for resolution of Acos# = 1/64. This figure also shows the 
tails local power index as a function of time (lower left) and 
the convergence order of the three resolution shown on the 
upper panel (lower right). 

In Fig. [3] we demonstrate the sixth order angular con- 



vergence of Kerr evolutions. This figure summarizes the 
main parts for our claim that our tail results are reliable: 
The code is convergent to a non-zero value, and even in- 
creases in magnitude with grid refinement (so that the 
signal seen is not an evolved noise signal), and the con- 
vergence order is global, throughout the computational 
domain. Notice that in Fig. [3] we study the convergence 
with the monopole projection of the total field. This 
projection is more sensitive to the created and evolved 
monopole numerical noise than the full field, say, as we 
can test its convergence also at times earlier than the in- 
stant the monopole dominates the total field for the first 
time. Figure [6] depicts the (L2 norm of the) full field 
and its monopole projection. The parameters are chosen 
so that intermediate tails are minimized, so that com- 
putation time is used efficiently. The monopole projec- 
tion amplitude is many order of magnitudes higher than 
typical "multipole leakage" noise levels (cf. with Fig. 
Indeed, it is the result of physical multipole excitations, 
not noise generation. Below, we also bring evidence that 
the observed tail is indeed asymptotic, so that we are 
not misled by an intermediate tail masquerading as an 
asymptotic one. 
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FIG. 3: Same as Fig. [1] for the monopole projection of the 
full field using the sixth order code for Kerr. The angular 
resolutions are Acos# = 1/80 (solid), 1/64 (dashed), and 
1/52 (dotted). The initial data are an outgoing Gaussian 
pulse located at r* = 25M with width of 4M with compact 
support on the Boyer-Lindquist t — hypersurface of a Kerr 
black hole with a/M = 0.995. The initial multipole is a pure 
Boyer-Lindquist £ = 4. The radial resolution is Ar = Af/50, 
and the temporal resolution is chosen so that At — Ar/2. 



III. THE TAIL IN AN a 2 EXPANSION 

GPP argue that when expanded in powers of a 2 , the 
tail would be at 0(a 4 ). When a mode P n exists at any 
iteration order (in a 2 ), at the next iteration three modes 
will result: P„_ 2 (if n - 2 > 0), P n , and P n+2 . The 
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FIG. 4: Same as Fig. [3] Shown are the L2 norm of the full 
field (solid curve) and the field's monopole projection (dashed 
curve). For both, the angular resolution is AcosS = 1/72 and 
the radial resolution is Ar = M/32. The temporal resolution 
is chosen so that At — Ar/2. 



mode excitation scheme is shown in Fig. [5] In the GPP 
approach, this behavior comes about because of the term 
sin 2 x P n (cos8) appearing in the source term (see [§| 
for detail). This way, starting with a pure £ = 4 mode, 
there is just one way to create a monopole field at 0(a 4 ), 
but three ways to create a monopole field at 0(a 6 ). We 
study two classes of initial data sets: In Class A the field 
ip is a pure multipole (in practice, we focus attention 
here mostly on the azimuthal £ — 4 mode) and the field's 
"momentum" II = 0, where II := dfip + bd? r ip = 0. 
Here, b = (f 2 + a 2 )/±, t 2 = (f 2 + a 2 ) 2 - a 2 A sin 2 6, 
and A = f 2 - 2Mf + a 2 (for details see, e.g., 0]). The 
f, 9 coordinates are the Boyer-Lindquist ones. Class A 
initial data lead to Tail A behavior (i.e., to n — 3), with 
the tail at 0(a 4 ). In Class B initial data we choose the 
field ip in the same way as in Class A, but replace the 
condition II = with dtip = ^df^ip or dtip = or 
a linear combination thereof (with coefficients that are 
independent of 9). Class B initial data lead to Tail B 
behavior (i.e., to n = 5), with the tail at 0(a A ). (The 
momentary stationary data case dtip — is unique, and 
leads to n — 6 in a similar way to the Schwarzschild case 
[To]].) The different sets of initial data are summarized in 
Table H 

From Fig. [5] is it evident why the late time tail for 
Class B data is at 0(a A ): The initial hexadecapole field 
at 0(a°) excites a quadrupole field at 0{a 2 ), which in 
itself excites a monopole field at 0(a 4 ). Higher powers 
in a 2 are also present in the monopole field, but for small 
values of a 2 they are negligibly small compared with the 
leading a 4 term. Even more so, the monopole tail appears 
to depend very weakly on terms at 0(a 6 ): the tail turns 



TABLE I: The initial data classes used in the text. Unac- 
cented coordinates mean any of the coordinate systems de- 
scribed in Table [TTJ See the text for more detail. 



Initial data class 


Condition 


Connection 


Class A 


at t = 0: ip, 11 = 


n = dtip + bdr,i> 


Class B 


at t = 0: ip, d t ip 





out to be a predominantly a 4 effect. 



P P P P 

2 4 6 



a 




a 

FIG. 5: The mode coupling schematics: at each iteration in 
a 2 in the GPP approach, each t modes excites three modes, 
specifically 1 — 2,1, and I + 2. The initial data are purely 
hexadecapolar. We show only modes up to 0(a 6 ) and £ = 6. 

Why then is the tail for Class A initial data also lead 
to a tail at (3(a 4 )? The condition n = is equiv- 
alent to the condition dtip = —bdf t ip. As the func- 
tion b depends on 9, it changes the multipolar struc- 
ture of the initial data sets. Specifically, for small a 2 , 
b= \ + {r-2M)a 2 sin 2 9/ (2f 3 ) + 0{a i ). One may there- 
fore naively expect these initial conditions to lead to a 
mixed state of multipoles at 0(a 2 ), that may include a 
monopole term at 0(a 2 ). If such a monopole term ex- 
isted, the tail would be an 0(a 2 ) effect. 

To find the order (in a 2 ) of the tail we write the 
monopole projection of the field ipo as an expansion in 
the spin parameter a: 

00 

flo = ^C„(i)xa". (1) 

71 = 

Class B initial data lead to only even powers of a having 
non-zero coefficients, and Co(t) — = C^it), so that the 
first non-zero coefficient is Ci(t). For the Tail B behavior 
found by GPP, Ci(t) ~ t -5 (for a pure Boyer-Lindquist 
initial mode with initial I = 4). 

For Class A initial data we can find the coefficients 
C n (t) from a 2+1D simulation in the following way: We 
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FIG. 6: The coefficients C 4 (£) and C 6 (t). The two sets of 
data points (o) corresponding to Ci(t) relate to fits to the 
ansatz Tp(t) — d,{t) a 4 (slightly higher) and ip(t) — C4(t) a 4 + 
Ce(t)a 6 . The coefficients Ce(t) are represented by □. We 
also show two reference lines, with slopes of —3 (solid) and of 
—5 (dashed). 



take the monopole projection value of the field in the 
tail regime at a fixed value of time for various a values. 
In practice, we considered a values in the range 0.05 < 
a/M < 0.995. We then fit the data to two ansatzs: first 
to -0o = C4 a 4 , and second, to tpo — C4 a 4 +Ca a 6 . That is 
we obtain the numerical values of the coefficients C n for 
that particular value of the time. The fit is very good 
we get the squared correlation coefficient to be better 
than 0.999. We then repeat for a number of time values, 
so that we have C±(t) and Cg(t). We can then find the 
time dependence of C±{t) and of C$(t). Our results are 
shown in Fig. [6l Both Ci{t) and C 6 (t) drop off as t~ 3 . 
We emphasize that the raw data are taken from Fig. [3] 
the data converge at sixth order angularly, and is the 
physical signal, not an evolved noise signal. 

We therefore make the following important inferences: 
first, initial data of Class A lead to to Tail A behavior. 
Second, the tail is a 0(a A ) phenomenon, and 0(a e ) con- 
tributions are negligibly small. We checked this result not 
only for small values of a 2 , but also for values as high as 
a 2 ~ 0.99. The reason for this behavior is as follows: the 
mixed state of multipoles at 0(a 2 ) in the initial data in- 
cludes only multipoles higher than the quadrupole. The 
absence of the monopole is evident from the vanishing of 



P (cos 0) x sin 2 P4 (cos ( 



sin 9d9 = 0. 



We do get however a quadrupole contribution at 0{a 2 ), 
as 



P 2 (cos0) x sin 2 0P 4 (cos < 



The monopole term on the initial data appears only at 
0(a 4 ). The function b is 



1 + 



f-2M 



a 4 sin 



2f 3 

2 a 



a 2 sin 2 1 



8f 6 

+ 0(a 6 ), 



3 cos 2 (f - 2Mf + (f - 2Mf 



16AP 



so that the term independent of inside the square brack- 
ets does not contribute to the monopole at 0(a 4 ), but the 
term proportional to cos 2 9 does, because 



P (cos 9) x sin 2 9 cos 2 9 P 4 (cos ( 



sin0d0 = 



16 

315 ' 



sin0d0 = -— ^ 0. 
105 



which does not vanish. 

We therefore find that the mixed state of multipoles is 
such that there is a monopole at 0(a 4 ), a quadrupole at 
0(a 2 ), and a hexadecapole at O(a ). These multipoles 
are arranged along a diagonal in Fig. [51 so that they 
track the mode coupling pattern, and therefore lead to a 
tail at 0(a 4 ), as each component of the mixed state of 
multipoles in itself leads to the same tail. Notably, when 
the initial data correspond to I — 6 the modes of the 
mixed state are again arranged along a diagonal, with 
a hexindatesserapole at O(a ), a hexadecapole at 0(a 2 ), 
a quadrupole at 0(a 4 ) and a monopole at 0(a 6 ). This 
behavior remains also for an initial £ = 8 data, with all 
mixed terms arranged along a diagonal in Fig. We 
therefore expect it to be a general result, and propose 
that Class A initial data tails are at the same order in a 2 
as Class B tails. 



IV. EVOLUTION OF CLASS B INITIAL DATA 
A. Axisymmetric modes 

Specifying initial data corresponding the a Class B hex- 
adecapole — namely, a pure £ = 4 moment on a Boyer- 
Lindquist time slice — we find the field and it's multipole 
projections to behave as in Fig. [7] Inspection of Fig. [7| 
shows we have not integrated long enough for the excited 
monopole to dominate the full field. We can neverthe- 
less still find the asymptotic tail index. This is done by 
finding the local power index of the monopole projection 
of the field (shown in Fig. [8]), and extrapolating it to in- 
finitely late times. We find the tail's slope to be 5.0129, 
within 0.3% from the value appropriate for Tail B. 

Specifying the Class B initial data to £ = 6, we fo- 
cus attention next on the subdominant quadrupole mode. 
Figure [5] depicts the full field and its hexadecapole and 
quadrupole moments, and Fig. 1101 shows the local power 
index for the quadrupole projection. Extrapolating to 
infinite times, we find the tail's slope to equal 9.0568, a 
deviation of 0.6% from the decay rate for the tail would 
be, if the quadrupole were the least multipole to be ex- 
cited. 
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FIG. 7: Same initial data as in Fig. [3] but for Class 
gaussian pulse is centered at r* = 25M with width 
Grid parameters are: AO = 1/48, Ar, = M/32. 
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FIG. 8: The local power index as a function of time for the 
monopole projection of the full field of Fig. [7] When extrap- 
olated to infinite time, the slope is = 5.0129. 
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FIG. 9: Evolution of Class B initial data for I = 6. Grid 
parameters are: A8 — 1/48, Ar = M/32. Shown are the 
full field (solid curve), and the hexadecapole (dashed curve) 
and quardupole (dotted curve) projections. (The monopole 
projection is not shown.) The initial data are again a gaussian 
centered at r* = 25M with width of 4M. 
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FIG. 10: The local power index as a function of time for 
the quadrupole projection of the full field of Fig. [9] When 
extrapolated to infinite time, the slope is n x = 9.0568. 



We therefore propose a generalization of the Tail B 
rule, that applies also for subdominant modes. Specif- 
ically, we propose that starting with an initial mode £, 
the decay rate of an (allowed) excited mode £' is given 
by n = I + £' + 1 if t - l> > 2, and n = 21' + 3 other- 
wise. This rule includes that of [7[ as a particular case, 
when £' = m, m + 1. As an example for the application 
of this generalized rule, consider the axisymmetric case 
above of £ = 6 and £' — 2 (which is a subdominant mode, 
as also the monopole is excited). The quadrupole mode 
£' = 2 then indeed decays with n = 9 as predicted by our 
proposal. 



B. Non— axisymmetric modes 

We next consider the case of non-axisymmetric initial 
data of Class B. Specifically, consider initial £ = 6 and the 
azimuthal mode is m = 2. In this case mode couplings 
allow the hexadecapole and quadrupole modes to be ex- 
cited, but the monopole excitement is disallowed. Ac- 
cording to our proposal, as £' = 2, £ — £' = 6 — 2 = 4>2, 
and therefore n = £ + £' + 1 = 6 + 2 + 1 = 9. We show in 
Fig. [TT] the full field and the quadrupole mode as func- 
tions of time, and in Fig. [T^] we plot the local power index 
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for the quadrupole projection as a function of (inverse) 
time. We therefore conclude that our proposal for the 
power law indices predicts the correct value also for non- 
axisymmctric initial data. To our knowledge, this is the 
first time that non-axisymmetric Kerr tails are consid- 
ered. 




FIG. 11: Evolution of Class B initial data for I = 6 and 
m = 2. Grid parameters are: AO = 1/48, Ar = M/32. 
Shown are the full field (solid curve), and the hexadecapole 
(dashed curve) and quardupole (dotted curve) projections. 
(The monopole projection is not shown.) The initial data are 
again a gaussian centered at r* = 25M with width of 4M. 
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FIG. 12: The local power index as a function of time for 
the quadrupole projection of the full field of Fig. [11] When 
extrapolated to infinite time, the slope IS 71 oo = 8.91. 
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FIG. 13: Solid curve: the full field. Solid reference line has a 
slope of 6.1. Dashed reference line has a slope of 4.9. Initial 
data were a Gaussian with 1 = 2. It was centered at 30. OM 
with a width of 4M. The grid resolution was: radial resolution 
Ar = M/20 and angular resolution Acos# = 1/32. Here, 
a/M = 0.8. 



V. OPTICAL ILLUSION 

When having asymptotic phenomena, it is crucial that 
the observations are made in the asymptotic regime. In 
Fig. [T3] we show the field starting at the same times, but 
extending to two different values of time. In the upper 
panel extending to t = 400M, and in the lower extending 
to t = 1600Af. In the upper panel we show a reference 
line with a slope of 6.1. One may be tempted to deduce 
from the seemingly two parallel curves at late times that 
the asymptotic slope of the field is 6.1. In the lower panel 
of Fig. [T3] we show the same field and the same reference 
line, but we add another reference line, with slope of 4.9. 
Extending the integration time, it now appears that the 
asymptotic slope is 4.9. Which is it then? 6.1 or 4.9? The 
correct answer is 3, as we specified here Class A initial 
data. The slope of 5.5 reported on in [f|, we believe, is 
just this effect: the field's decay rate was determined in 
5] not in its true asymptotic regime. 

In Fig. [U we show the full field up to t = 500Af and 
to t/M = 2950. At first we show a reference curve with 
slope of 5.5, and then increase the integration time, and 
show both the first reference curve, and a new one with 
slope of 3.0. In particular, if we were uncareful we could 
report on a tail index of 5.5 like that of Continuing 
the simulation, however, we find the true asymptotic tail, 
with a tail index of 3, as expected. 

The arbitrariness of the decay rate, as determined 
from intermediate time segments, is further exemplified 
in Fig. [151 which shows the field for a number of dif- 
ferent choices of tq, the center of the initial pulse. The 
closer the pulse is to the black hole, the sooner the true 
asymptotic tail appears. Conversely, pushing the initial 
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FIG. 14: Solid curve: the full field. Solid reference line has a 
slope of 5.5. Dashed reference line has slope of 3. Initial data 
were a Gaussian with 1 = 2. It was centered at 17. 5M with 
a width of 4M. The grid resolution was: radial resolution 
Ar = M/20 and angular resolution A cos (9 = 1/32. Here, 
a/M = 0.8. 




FIG. 15: Effect of varying the location of the initial pulse. In 
all these runs, initial data was a Gaussian of width A = 4Af, 
with 1 = 2, and a/M = 0.8. The grid resolution was: 0.05 (r) 
x 0.03125 (angle). The curves shown and the corresponding 
value for ro, the center of the initial data gaussian, are: dotted 
curve: ro = 10M, dash-dotted curve: ro = 12.5A/, dashed 
curve: ro = 15M, solid curve: ro = 17. 5M, dashed thick 
curve: ro = 20Af, and solid thick curve: ro = 30M. All data 
are extracted at r = 17. 5M. 

pulse outwards, we protract the intermediate tail regime, 
and can make it as long as we want: it is fully control- 
lable. The choice of the length of the intermediate tail 
also affects its decay rate determination. We therefore 
see that not only is such a decay rate meaningless, it is 
also directly manipulable. The important lesson here, re- 



flecting on the claims made in [5j , is that a distant initial 
pulse makes it harder to see the asymptotic tail, and it 
would require longer time evolutions. In the initial 
pulse was chosen at ro = 100M, which would make the 
intermediate tail very protracted indeed. 




FIG. 16: Same data as in Fig. 1131 The full field is decomposed 
into the Boyer-Lindquist multipole modes. 

Figure [TB] explains the broken power-law behavior we 
see in Fig. 1131 and also why neither is the true asymp- 
totic field: Up to t ~ 500M the field is dominated by 
the dipole (£ — 2) mode, and at later times it is dom- 
inated by the monopole (£ = 0). Notice, that even the 
seemingly constant slope in the lower panel of Fig. [13] 
does not represent the asymptotic behavior of the field, 
as is clear from Fig. [TH] The "optical illusion" effect is 
not a result of mode mixing: it appears in each mode 
separately. Indeed, in the t = mode what may appear 
as a straight segment (on the log-log plot of Fig. [T6|) is 
curving if the integration time is extended. But when we 
combine the two separate optical illusions of the 1=1 
and i = modes, we may get a broken power-law as in 
the lower panel of Fig. [131 that has nothing to do with 
the true asymptotic decay rate of the field. Notice, that 
also the i = 4 is excited, but at a negligible level. Higher 
(even) modes are also excited, but not shown in Fig. 1161 

In view of he arbitrariness in the slope as calculated 
over some finite range, it appears that a better way to 
evaluate the slope is through its local value, the local 
power index defined by n(t) :— —ipt/iji Even plotting 
the local power index as a function of the time could be 
misleading, as is evident from Figs. [T7] and |TS] notice 
how the bold solid curve appears to approach a constant 
value for low values of t. Indeed, changing the initial 
conditions, one can made the local power index have an 
arbitrarily protracted intermediate value. 

How can then one distinguish an arbitrarily protracted 
intermediate behavior form the true asymptotic behav- 
ior? The key point in distinguishing the true asymptotic 
tail from an intermediate one is to consider the asymp- 
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data pulse, so that one has a protracted intermediate tail, 
and the other arrives quickly at the asymptotic regime. 
Indeed, not only does the squared correlation coefficient 
indicate only relatively poor agreement with the ansatz 
([2]) for the intermediate tail, it also fails to improve at 
later times. In contrast, at the true asymptotic regime, 
the squared correlation coefficient approaches unity con- 
sistently. We are therefore motivated to introduce the 
following definition: 

Definition: The local power index satisfies the ansatz 
@ locally if the squared correlation coefficient for fitting 
the local power indices of constant duration intervals to 
the ansatz ([2]) approaches unity as t increases. 
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FIG. 17: Local power index (LPI) as a function of t for the 
same data as in Fig. [15] for the I — Boyer-Lindquist multi- 
pole mode. 
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FIG. 19: Dependence of the squared correlation coefficient R 
on the local temporal location of an interval of constant lapse 
(lasting lOOOAf). Data are taken for the same initial data as 
in Fig. 1151 For both cases the data correspond to the squared 
correlation coefficient of ten evenly spaced (in t/M) values of 
the local power index (with increments of 100A/) ending at the 
value shown, fitted to the ansatz ([2}: Squares (□) correspond 
to data for ro = 30M, and circles (o) to ro = 10M. 



FIG. 18: Same as Fig. 1171 but presented as a function of M/t. 

totic behavior of the local power index. As shown in [ll[ , 
the asymptotic behavior of the local power index is 

n(t) =n QO +n 1 (M/t) + ■■■ . (2) 

One should therefore test how well the fit of n(t) as a 
function of M/t improves with t. Loosely speaking, in- 
termediate regions, like those seen in Fig. 1181 will gener- 
ally yield fits that deteriorate with time, until the true 
asymptotic regime is reached. How should this criterion 
be applied in practice? A necessary condition is that the 
ansatz is satisfied locally. Specifically, Fig. [TH dis- 
plays the squared correlation coefficient 1 — R 2 for con- 
stant duration intervals (of ten equally spaced in time 
values of the local power index). In practice, we choose 
two sets of data differing only in the location of the initial 



The local satisfaction of the ansatz ^ is only a nec- 
essary condition for the true asymptotic tail, because in 
practice we can only examine numerically finite time val- 
ues. Indeed, this is the very reason why misidentifying 
an intermediate tail as an asymptotic one is a problem 
in the first place. Perhaps not surprisingly, one can find 
intermediate tails for which the ansatz |(5J) appears to be 
satisfied locally, but is violated globally in the following 
sense. 

Definition: The local power index satisfies the ansatz 
([2]) globally if the squared correlation coefficient for fitting 
the local power indices of consecutively longer intervals 
to the ansatz @ starting from an arbitrary (late) time 
approaches unity as t increases. 

In Fig. [20] we show the local satisfaction but global 
violation of ansatz ^ for an intermediate tail (upper 
panel ([2"DR )). and their simultaneous satisfaction for the 
asymptotic tail (lower panel (|20B)). 
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We are therefore motivated to propose the following as 
a criterion for distinguishing intermediate from asymp- 
totic tails: 

A criterion for identification of the asymptotic 
tail: The tail is asymptotic if the ansatz ([2]) is satisfied 
simultaneously both locally and globally. 
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FIG. 20: Dependence of the squared correlation coefficient R 2 
on a global temporal interval. Upper panel (A): behavior of 
an intermediate tail. Lower panel (B): behavior of an asymp- 
totic tail. In both panels the local satisfaction of ansatz © 
is shown with circles (o), and the global violation (A) or sat- 
isfaction (B) of the same ansatz is shown with squares (□). 
Initial data are centered on ro = 30Af (same as in Fig. I15[l . 
Data for the upper panel (A) correspond to five evenly spaced 
(in t/M) values of the local power index (with increments of 
30M) ending at the value shown (local) and appropriately 
many values with the same spacing (starting at t — 360A/) 
and ending at the shown time. Data for the lower panel (B) 
correspond to five evenly spaced (in t/M) values of the local 
power index (with increments of 100M) ending at the value 
shown (local) and appropriately many values with the same 
spacing (starting at t — 500A/) and ending at the shown time. 

While figures such as Fig. 02] or the upper panel of 
Fig. [T3]may lead to a misidentification of an intermediate 
tail as an asymptotic one, we believe that our criterion 
will expose an intermediate tail as such, and correctly 
identify asymptotic tails. 



VI. INGOING KERR SLICING 

The tails for initial data specified on ingoing Kerr slices 
were first discussed in [3| • The initial data used in Q are 
of Class A, leading to Tail A behavior. Recently, it has 
argued that a pure multipole in the ingoing Kerr coordi- 
nates does not correspond to a pure multipole in Boyer- 
Lindquist coordinates, and that in addition the different 
slicing conditions in the two cases were the reason that 
[|| found Tail A behavior and not Tail B behavior. 



Ingoing Kerr coordinates are different from Boyer- 
Lindquist coordinates by a transformation of the <fr and 
t coordinates — with the transformation depending on 
the radial coordinate only — leaving the ingoing Kerr 
r, coordinates identical to their Boyer-Lindquist coun- 
terparts. Specifically, the transformation is 4> = 4> + 
a J A -1 df and t = t + r* — f, where t, <fi are the ingo- 
ing Kerr coordinates, t, <f> are the Boyer-Lindquist co- 
ordinates, and where A = f 2 — 2Mf + a 2 and r* = 
j(f 2 + a 2 )A _1 (if. The different sets of coordinates we 
use are summarized in Tabic ILTl 

TABLE II: The notation used for the different coordinate sys- 
tems and the transformations between them. See the text for 
more detail. 



Coordinates 


Notation 


Transformation 


Boyer-Lindquist 


i, f, 9, 4> 




Ingoing Kerr 


t,r,9,<j> 


= + a J A _1 df 
t = t + r* — f 


Poisson 


t, f, 9, 4> 


f = \] r 1 + e 2 sin 2 9 

r^-j-e sin* tf 



A spherical harmonic mode in Boyer-Lindquist coordi- 
nates with a specific multipole £ and a specific azimuthal 
mode m is therefore also a spherical harmonic mode in 
ingoing Kerr coordinates with the same values for £, m, 
and the difference is only in a different function of the 
radial coordinate. One may therefore expect that initial 
data specified in ingoing Kerr coordinates would lead to 
the same tail behavior as in Boyer-Lindquist coordinates. 
One may object to this expectation by stating that the 
constant t slices in the Boyer-Lindquist case are not the 
same as the constant t slices in the ingoing Kerr case. 
While this difference is certainly existing and the two 
hypersurfaces are not the same, we find that it does not 
change the tail behavior. 

Class A initial data were studied in detail in We 
therefore present here the results for Class B initial data. 
We find that starting with a pure I — 4 azimuthal 
mode (in ingoing Kerr coordinates onat = const slice) 
the tail has Tail B behavior, similar to that of Boyer- 
Lindquist coordinates and Boyer-Lindquist slicing. Fig- 
ure shows the field and the local power index for these 
data. We infer from the simulation data that the tail 
index is 5.088, a deviation of 1.7% from the Tail B pre- 
diction of 5. 

The ingoing Kerr results are obtained using a code 
based on the Penetrating Teukolsky Code (PTC) [lj]. 
The PTC is written on a fixed ingoing Kerr coordinate 
grid, which means that grid resolution becomes increas- 
ingly poor approaching the event horizon. While simple 
Fixed Mesh Refinement methods can resolve this diffi- 
culty, we were able to obtain the tail index using the 
following method: we first find the local power index for 
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FIG. 21: Evolution of Class B initial data of a pure azimuthal 
scalar field hexadecapole (£ — 4) mode in ingoing Kerr coor- 
dinates on a t = const slice. Here, a/M — 0.8. The initial 
data radial profile is a gaussian with ro = 2M and A = 1M. 
The grid resolution is A8 = 1/64 and Ar = 2M/125. The 
dashed curve is the carve ip = (M/t) which is added for ref- 
erence. The local power index is calculated for three radial 
resolutions: Ar = M/40 (dotted), Ar = M/50 (dashed), and 
Ar = 2M/125 (solid). 



a number of time values for different choices of the ra- 
dial (and temporal) resolution. We then extrapolate the 
local power indices to infinite radial resolution, and then 
extrapolate the sequence of local power indices at dif- 
ferent values of the time (and with infinite resolutions) 
to infinitely late times. The result is rather close to the 
expected value of 5, but suffers from somewhat higher 
numerical error because of the two extrapolation cycles. 
As noted above, the deviation of the tail index from the 
Tail B value is at the order of 10 -2 . We therefore con- 
clude that Class B initial data in the ingoing Kerr case 
lead to Tail B behavior. 

We therefore contend that initial time slices, even when 
are distinct, fall into equivalency classes in terms of 
the type of tail behavior that emerges (recall that on 
a given slice multipolcs can also mix because of coordi- 
nate transformations). The Boyer— Lindquist and ingoing 
Kerr slices belong to the same equivalency class, in the 
sense that one finds the same Tail B behavior for Class 
B initial data. The Kerr-Schild slices belong to a dis- 
tinct equivalency class, and Class B initial data lead to 
Tail A behavior. We conjecture that the Boyer-Lindquist 
and ingoing Kerr slices belong to the same equivalency 
class because the time slices in either are related by a 
transformation that does not involve the polar angular 
coordinate and that depends only on the temporal and 
azimuthal angular coordinate. One may therefore con- 
jecture that other slicing conditions that differ from the 
Boyer-Lindquist slicing condition by similar transforma- 



tions will behave accordingly. 



VII. POISSON DATA 

Poisson has made the argument that if Class B initial 
data are specified in so-called "spherical coordinates," 
rather than in Boyer-Lindquist coordinates, which Pois- 
son calls "spheroidal," mode evolution is changed [6J. In 
particular, Poisson proposes that in Kerr spacctime such 
initial data lead to Tail B behavior. In a spacetime with 
global weak curvature (no black hole, but no assumptions 
on symmetry), Poisson data indeed lead to such late- 
time behavior. This may be surprising, as a pure mode 
of Poisson coordinates may be written as a mixed state 
of Boyer-Lindquist modes. This mixed state of modes 
effectively behaves as Class A initial data, so that the 
intuitively expected behavior is that of Tail A. 

Clearly, the assumption of globally weak curvature af- 
fects the tails, and shifts them from their Tail A values to 
those reported by Poisson. This could be the case if the 
field reflects off the regular center in the Poisson space- 
time, so that the leading order part of the field (in i -1 ) is 
cancelled out and one is left with the next order in t _1 , 
or if the monopole part of the field is canceled after the 
reflection, leaving only the dipole and higher moments 
present Q. The following question arises: what are the 
tail indices for "pure Poisson" Class B data in a Kerr 
spacetime? Specifically, if one sets at t = a pure Class 
B hexadecapole mode with respect to the Poisson "spher- 
ical" coordinates, does the late-time tail decay as 1/i 3 , 
as appropriate for Tail A, or as 1/t 5 as appropriate for 
Tail B (and as predicted by @)? As we are showing here, 
the answer is clearly 1/t 3 . 

Initial data that describe a "pure Poisson" mode are 
a linear combination of Boyer-Lindquist modes, as one 
may expand such a mode in the function space of Boyer- 
Lindquist spherical harmonics. Importantly, as the trans- 
formation from Boyer-Lindquist to Poisson coordinates 
does not involve the time, Boyer-Lindquist slices are 
identical to Poisson slices. Therefore, one may use a code 
in which the wave equation is written in Boyer-Lindquist 
coordinates to simulate the evolution of a "pure Poisson" 
mode, if one writes the initial data as the appropriate lin- 
ear combination of Boyer-Lindquist modes. This may be 
done as follows: Denote the Boyer-Lindquist coordinates 
as f, 6 and the Poisson "spherical" coordinates as r, 6. As 
shown in |h|, the coordinate transformation is 



cos 2 fl = 



r 2 + e 2 sin 2 i 



f 2 + e 2 sin 2 I 



(3) 
(4) 



where e is a measure of the "ellipticity" of the Boyer- 
Lindquist coordinates, which in practice equals the Kerr 
black hole's spin a. To specify a "pure Poisson" mode, 
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FIG. 22: A scalar field "pure Poisson" i = 4 mode, C^ff) 
(solid curve), and its projections on Boyer-Lindquist £ modes 
4 CV(f): dashed (£' = 4), dash-dotted (£' = 2), and dotted 
(£' = 0). Here, a/M = 0.995. The initial data radial profile is 
a gaussian with fo — 25M and A = 4M. The abscissa is r/M 
or f/M as is appropriate. 



A, with the monopole's decay rate at infinite time ex- 
trapolated to 3.0017, within 6 x 10 -4 from the predicted 
value of 3. 



o 




iog 10 f 
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FIG. 23: Time evolution of a Class B "pure Poisson" £ = 4 
mode, and its Boyer-Lindquist £' — 0, 2, 4 modes, and the lo- 
cal power indices for the full field and for the Boyer-Lindquist 
monopole projection. Here, a/M = 0.995. The (Class B) ini- 
tial data are those of Fig. [22] 



we take the field to be 

1>t\t=Q = C t (?) P t (.COB§) 

= Acxp ^^ 2 ) P^cosfl), (5) 

after choosing the radial function to be a gaussian of 
width A centered on fo. The initial "pure Poisson" field 
V^lt=n can be expanded in the function space of Boyer- 
Lindquist Legendre polynomials as 

oo 

M=o = £ *-df{f)Pf {cos6). (6) 

l'=0 

As shown in the functions iCv(f) are complicated. 
(In fact, only the coefficients for the inverse transforma- 
tion are found in [6( (for £ = 4 and £' = 0) as a power 
series in the "ellipticity" e.) In Fig. [22] we show the func- 
tions C±(r), and iCi'{f) for £' = 0,2,4 (in general, these 
functions are non- vanishing for all values of £' ', respect- 
ing the parity symmetry), for a high value for the "el- 
lipticity," specifically for e = 0.995A/ . To complete the 
initial data problem, we also need to specify at t = 
the time derivative of the field. For Class B initial data 
we choose dtip = df,ip, and for Class A initial data we 
choose II — 0. 

Figure[23ldepicts the full field and its multipole projec- 
tions and the local power indices for the full field and its 
monopole projection as functions of time for Class B ini- 
tial data. The late time behavior is clearly that of Tail 



In Fig. [24] we consider Class A initial data. Again, it 
suggests that the tail field drops of as l/t 3 for a scalar 
field with initial data that are a "pure Poisson" hexade- 
capole mode. Specifically, we show this to be the case 
more precisely, using the criterion based on simultaneous 
local and global behavior of the local power index. Com- 
parison of the evolution of the evolution of the "pure 
Poisson" £ = 4 mode with the evolution of its Boyer- 
Lindquist £' = 4 projection suggests that the monopole 
projection £' = at the initial time is important at late 
times. The £' = mode is being excited during the evo- 
lution also when it is not present at the initial time, and 
eventually it dominates the field at late times. However, 
the generated £' = mode starts dominating the total 
field only at a later time, as its amplitude is lower than 
that of the monopole part of the full £ = 4 mode. The 
latter mode includes a monopole component already at 
the initial time, and this part of the field starts dominat- 
ing the total field at an earlier time than the monopole 
component that is excited. For this reason the asymp- 
totic £ = 4 mode is much stronger than the asymptotic 
tail of initial data that include only the £' = 4 projection 
of the £ = 4 mode, even though at the initial time they 
are almost identical, as can be seen in Fig. [22] 

To show that indeed the late-time tail is dominated 
by the evolution of the £' = projection of the original 
"pure Poisson" £ = 4 mode at t — 0, we present in Fig. [231 
the evolution of this mode. Comparison with Fig. [24] in- 
deed shows that the late time field of the full £ = 4 mode 
is dominated by the evolution of its monopole projec- 
tion, not the monopole mode that is excited during the 
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FIG. 24: Time evolution of a Class A "pure Poisson" £ — 4 
mode (panel A), and its Boyer-Lindquist £' = 0,2,4 modes 
(panels B,C, and D, respectively) (solid curves), and the evo- 
lution of the Boyer-Lindquist £' = 4 projection of the "pure 
Poisson" mode and its Boyer-Lindquist projections (dashed 
curves). Here, a/M = 0.995. The initial data are those of 
Fig. m 
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FIG. 25: Time evolution of the Boyer-Lindquist £' — pro- 
jection of the Class A "pure Poisson" t = 4 mode. The initial 
data are those of Fig. I22I 



FIG. 26: The local power index as a function of t~ x for the 
Class A "pure Poisson" mode I = 4 (solid curve) and for its 
£' = Boyer-Lindquist projection (dashed curve). The initial 
data are those of Fig. I221 and the field evolution is presented 
in Fig. [21 



the initial data of a Y^^° mode, the i^™^ mode is 
excited. Figure [24] (in conjunction with the conclusion 
that the tail is asymptotic) shows that when starting 
with a "pure Poisson" mode, the excited monopole mode 
is insignificant for the asymptotic tail problem, as it 
is swamped by the evolution of the monopole projec- 
tion of the initial data. Above, we showed that start- 
ing with Class B initial data of a pure Boyer-Lindquist 
I = 4, to = 0, the late time tail drops off according to t~ 5 , 
in accordance with Tail B behavior, consistently with the 
results of [§, [l(| • In that case, there is no monopole field 
present at t = 0, and all the monopole field present at 
late times is due to the excitation of the £ — mode due 
to mode coupling. The excited mode that has Tail B be- 
havior is overwhelmed by the evolution of the monopole 
mode existing at the initial time that has Tail A behav- 
ior. Most importantly, we choose the initial data to be 
a pure spherical harmonic mode in the so-called "spheri- 
cal" coordinates. The late time tail still decays according 
to Tail A behavior, in contrast with the Tail B expecta- 
tions based on [fj. In the case presented, as the initial 
mode is azimuthal, the asymptotic tail drops off as t~ 3 . 



evolution because of mode coupling. 

In Fig. [26] we show the behavior of the local power 
indices for the full field of the Class A "pure Poisson" 
£ = 4 mode, and for its monopole £' — projection. 
Figure [26] strongly suggests that the full field, and also 
the monopole projection, drop off at late times as 1/t 3 . 
Indeed, the tail seen in Figs. [25] and [26] is asymptotic, 
and the local power index satisfies the ansatz both 
locally and globally simultaneously. 

As noted in [||, it is not surprising that starting with 
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